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Abstract — Aggressive frequency reuse in the return link (RL) 
of multibeam satellite communications (SatComs) is crucial to- 
wards the implementation of next generation, interactive satellite 
services. In this direction, multiuser detection has shown great 
potential in mitigating the increased intrasystem interferences, 
induced by a tight spectrum reuse. Herein we present an analytic 
framework to describe the linear Minimum Mean Square Error 
(MMSE) performance of multiuser channels that exhibit full 
receive correlation: an inherent attribute of the RL of multibeam 
SatComs. Analytic, tight approximations on the MMSE perfor- 
mance are proposed for cases where closed form solutions are not 
available in the existing literature. The proposed framework is 
generic, thus providing a generalized solution straightforwardly 
extendable to various fading models over channels that exhibit 
full receive correlation. Simulation results are provided to show 
the tightness of the proposed approximation with respect to the 
available transmit power. 

Index Terms — Multiuser detection, multibeam satellites, return 
link, linear minimum mean square error receivers. 



I. Introduction & Related Work 

Towards next generation broadband, interactive, multibeam 
satellite systems, higher throughput requirements especially 
for the return link (RL) are necessary. Multi user detection 
(MUD) has already proven a key technique for the enhance- 
ment of the spectral efficiency of the the RL (£TJ and the 
Forward Link (FL) [2 1 of multibeam satellite communications 
(SatComs). However, the inherent differences between the 
satellite and the terrestrial channel impose difficulties in the 
theoretical performance analysis of MUD for satellite systems. 
Added to that, linear MMSE receivers are more difficult to 
analyze compared to the optimal capacity achieving successive 
interference cancellation (SIC) techniques. 

The basic characteristic of the RL multibeam satellite chan- 
nel is the high receive correlation among the channels at the 
satellite side, resulting from the total lack of scatterers, as 
well as the practical collocation of the multibeam antenna 
feeds, at the satellite side. As a result, all the receive antennas 
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experience identical fading instances from a specific user. 
Specifically, each user's vector channel towards the satellite 
antennas is defined by the antenna radiation pattern. According 
to the multiplicative model and assuming block fading, the 
product between a single fading coefficient and the channel 
vector will characterize the fading instance. Subsequently, the 
channel matrix is given by the product of a full rank multibeam 
gain matrix with a diagonal matrix, composed of random 
fading coefficients, as originally proposed in [1| and extended 
in [3 1 and [4|, for more realistic system models. 

Linear multiuser receivers where extensively examined in 
0. More recently, numerous contributions examined the 
asymptotic performance of these architectures by employing 
random matrix theory methods described in [6| and the refer- 
ences therein. In the same field, the work of Q is also noted. 

Moreover, the mutual information of channels for which 
exact analytical formulas do not exist in literature, has 
been described via analytical, tight bounds. To this end, the 
Minkowski's inequality was initially proposed in [8 1 to provide 
a tight lower bound on the channel mutual information. 
This technique was extended in [9| to provide a new lower 
bound on the ergodic capacity of distributed MIMO systems. 
Furthermore, generic bounds on the sum rate of zero forcing 
(ZF) receivers over Rayleigh, Rice and correlated channels 
were given in [ 10 1. 

Focusing on SatComs, the performance of linear MMSE 
receivers is not trivial to characterize analytically. In the 
existing literature, two different modeling approaches have 
been considered. The authors of ifTTI have modeled fading 
as the sum of two random matrices, each following a specific 
distribution, without however reporting results for linear re- 
ceivers. The other approach, namely the multiplicative model, 
was introduced in JT], where the performance of optimal 
non-linear and linear receivers was given as a function of 
matrix arguments, with high computational complexity, for 
theoretical channel models. This model has been extended in 
for composite, realistic channels and a novel, simple lower 
bound for the channel capacity has been deduced. Finally, the 
incorporation of rain fading has been tackled in [4|, using 
an identical channel model, in which asymptotic closed form 
expressions for the system's capacity and the MMSE were 
derived. 

The purpose of the present contribution is to extend the 
results of J3] and E) for the case of linear MMSE receivers, 
thus providing an analytic framework that covers the gen- 
eralized type of channels with full receive correlation. To 
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this end, the information theoretic link between the channel 
mutual information and the MMSE given in [ 12 1, is exploitecQ. 
Results on the MMSE performance of the generalized satellite 
channel, for which the mutual information has been analyti- 
cally described only through tight bounds, are proposed. 

The rest of this paper is structured as follows. In Sec. [TT] 
the generalized multibeam satellite channel is introduced. The 
analytic framework proposed is described in detail in Sec. 
Hill Simulation results are presented in Sec. IIVI while Sec. 
concludes the findings of this paper. 

Notation: Throughout the paper, £{■}, (■) t , denote the 
expectation, and the conjugate transpose of a matrix. Bold 
face lower case characters denote column vectors and upper 
case denote matrices, while I„ denotes an identity matrix of 
size n. 

II. Signal & Channel Model 

Let us consider a multi-user (MU) single input multi- 
ple output (SIMO) multiple access channel (MAC) with K 
single-antenna terminals transmitting towards a single receiver 
equipped with N = K antennas^. The input-output relation- 
ship reads as 

y = ^Hs + n, (1) 

where s 6 C is the transmitted baseband signal vector, 
such that £{ss ff } = Ik, y € C Nxl is the received signal 
vector, n ~ CAf (0,1) is the complex noise vector and 7 
is the normalized transmit signal-to-noise ratio (SNR). Also, 
the matrix H 6 q NxK represents the complex- valued block 
fading channel and admits a number of expressions depending 
on the model. 

A. Generalized Multibeam Channel: In the present contri- 
bution, the channel matrix is assumed to be modeled as the 
product of a full rank matrix with a random fading matrix. The 
assumption of full receive correlation, that will be hereafter 
considered, shall impose a diagonal fading matrix 0], on 
the multiplicative model, thus yielding 



H = BD 1/2 . 



(2) 



where B is a fixed, deterministic, full rank matrix with 
normalized power (i.e. tracc(B 2 ) = A") and D = diagQ is a 
random diagonal matrix composed of the fading coefficients. 
The deterministic matrix B models the multibeam antenna 
radiation pattern, while the random elements of the diagonal 
matrix D, are drawn from distributions that model various 
small or large scale effects. Subsequently, the adopted model, 
provides a generalized approach towards systems with full 
receive side correlation. By changing the fading distribution 
different systems can be modeled. Two specific examples, 
relevant for multibeam SatComs will be studied hereafter. 

1 For completeness we also denote the work of | 1 3 | that provides an analytic 
framework for linear MMSE receivers, which however does not cover our 
case. 

2 The commonly adopted assumption of symmetric systems (e.g. (9)) is 
a requisite for the analysis. In multibeam SatComs, it models systems that 
schedule a single user per beam, each time slot. This assumption is usually 
adopted in the relevant literature (T), (5J, 1141 . 



1 ) Composite fading: This channel has been proposed in 
to model the multibeam mobile satellite channel, by incor- 
porating small scale Rician fading and large scale lognormal 
shadowing. The resulting model reads as 



H r — BHfj 



(3) 



where ILj and X<j are diagonal fading matrices with elements 
drawn from Rice and log-normal distributions respectively. 

2) Rain fading: When fixed satellite services are assumed, 
then the high antenna directivity imposes an AWGN channel 
and allows for the utilization of higher frequency bands, where 
rain fading can dramatically deteriorate system performance. 
For this case the model proposed in [4| will yield 



H, 



(4) 



where the random elements of d are drawn from a log-normal 
distribution in dB scal^. 



III. Analysis 

The performance of linear multiuser detection (MUD) is 
evaluated via the achievable SINR^ after MMSE detection at 
the fcth user, given by |[T6ll 
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/mmse,& 



(Ijc + 7HtH)- 



- 1. 



(5) 
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Averaged over the users and all channel instances, the system 
spectral efficiency, will be given by lfl6l 



C mmse =£h{^E l0g 2 ( X + 7mmse, k) j , 



(6) 



in bits/sec/Hz. By combining (O and (|6]l and directly applying 
the Jensen's inequality the following stands 



nse > £h j - log a (j£ trace ((l K + yftfH) 



(7) 

Thus, in common practice, another measure, namely the in- 
stantaneous, per user, minimum mean square error, is also 
adopted l6l. fl2l: 



i trace (jl K + 7 H f H) 1 



(8) 



The presence of the inverse of a matrix sum in (jTJ hinders 
the computation of e 2 when the eigenvalue distribution of 
H^H cannot be computed analytically. To solve this problem, 
we propose to use the approximation that follows. 

Lemma 1. The per user MMSE of an uplink MU SIMO system 
can be approximated by 

1 



l + 7exp(^lndet(HtH))' 
relative to a specific channel instance. 



(9) 



3 In natural units, these elements are produced by exponentiated log-normal 
elements 1151 and the relevant distribution is also referred to as log-log- 
normal. 



Proof: In fl2l . an explicit relationship between the 
channel mutual information I e = log 2 det (I + 7H^H) in 
bits/sec/Hz and e 2 was derived by differentiating with respect 
to the SNR. This result is easily extended for vector channels 
yielding [6| 



7^-I e (7) = K- trace + 7H t H) 



K(l-e 2 ). 



(10) 
(11) 



Furthermore, with the use of Minkowski's inequality, a bound 
on I e , tight over the whole SNR region and exact for the high 
SNRs, has been derived in (S): 



l e >hb =Klog 2 (l + 7 det((HtH 



.l/K 



(12) 



By differentiation with respect to 7 we get the following: 



7 Xdet((HtH) 1/A ' 



1—hb = 
9 t" l + 7 det((HtH) 1/K 



(10) 



e 2 = 



1 



1 + 7 det ( (HtH 



,1/K 



(13) 



Finally ( fT3l can be rewritten as in ([9j. ■ 
Since the differentiation does not preserve the direction of 
the bound, the characteristics of this approximation will be 
studied in more detail in the following lemma. 

Lemma 2. Let e 2 be the derived approximation of e 2 , then e 2 
and e 2 as functions of 7, will have a crossing point. 

Proof: Denoting as Aj the ordered eigenvalues of 
H^H, let us define the function D (7) = I e (7) — 
I«,(7), where J e (7) = £\ log (1 + 7^) and I lb (7) = 
log ^1 + 7 (Yii -^i) 1 ^) ' w i m the following properties: 1) 
-D(O) = 0. 2) For -f a sufficiently large, but finite, we have 
that D(j ) = 0. From properties 1 and 2, a straightforward 
application of Cauchy's mean- value theorem yields that there 
will be at least one point 7* 6 (0, 7 ), with zero derivative 
and subsequently a crossing point between the approximation 
and the actual function. ■ 

Lemma 3. The average over the channel realizations, per user 
MMSE, £ H {e 2 }, can be bounded by 



(14) 



<cxp (±£ H I In det (HtH) 



1 



Proof: Let us consider the function <p(x) = 
((1 + exp (x))) 1 which is convex for x < and concave for 
x > 0. By applying Jensen's inequality over the two regions 
we get that 

£h{e 2 } = <! n (15) 

I < a, x > 

for a = (\ + 7 exp j In det (H^H) and 

^ = l/K-lndet((HtH))+lri7. ' ■ 





Figure 1 . Multibeam antenna beam pattern for European coverage designed 
by the ESA over which a specific cluster of 7 beams (in red), is chosen. 



Lemma 4. The average per user MMSE approximation, for 
the composite multibeam satellite channel (as given by (fjjj, is 
analytically described by ( 1181 l, where fx m (dBj is the mean of 
the normal distribution, K r the Rician factor, the §2 function 
is given as <?2(s 2 ) = logs 2 + Ei(s 2 ) (as defined in /E1?J, Ei 
is the exponential integral and s 2 — K r is the non-centrality 
parameter of the associated x 2 -distribution. 

Proof: The expectation of the logarithm of the determi- 
nant of the specific channel matrix has been derived in Q. 
Direct application of these calculations on (fl4l concludes the 
proof. ■ 

Lemma 5. The per user average MMSE spectral efficiency of 
a multibeam satellite system under rain fading (as modeled by 
can be approximated by the closed form formula 

1 



1 + 7 exp ( -iln(dctB 2 ) 



(16) 



where /ii the mean of the equivalent log-normal distribution. 

Proof: Following the same line of reasoning as in the 
proof of Lemma 3, an analytical bound on £n {e 2 ain } will 
read as 

(17) 

l+7exp(i£ H {lndet(B 2 d )}) 

where lndct(B 2 c ;) = In (detB 2 ) +lndetd since the matrices 

are square and £ {In (det d)} = X^i m ^ = K ' Mi- The 
parameter p,\ is the mean of the related log-normal distribution 
given by fn — exp(/i m + c 2 /2). As before, fx m and a are the 
mean and variance of the associated normal distribution. Since 
B is a deterministic matrix, from the above analysis, dT6l is 
deduced. ■ 

IV. Simulation Results 

Realistic beam patterns from multibeam coverage over 
Europe, that is a 245-beam pattern designed by the Euro- 
pean Space Agency (ESA), are chosen to provide accurate 
results for the system performance (Fig. Q}. In full-frequency 
multibeam systems, multiple GWs are deployed due to limited 
feeder link bandwidth. Each GW serves only a small number 
of beams. Therefore, only a part of the total beam pattern (i.e. 
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1 + 7 exp i log (det B 2 ) + t i m + g 2 (s 2 ) - log (K r + 
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Figure 2. MMSE for 7 normalized ESA beams with K r = lOdB, /i„ 
-2.63, a=0.5. 



7 beams) has been considered in the analysis as presented 
in Fig. Q] These beams correspond to the beams served by a 
specific gateway (GW) over which MUD will be performed 
in a distributed GW scenario. Monte Carlo (MC) simulations 
have been performed to calculate the exact performance over 
the random channel. The analytic results of Sec. [Ill] provide 
deterministic approximations on the system performance with 
respect to 7 and are compared to the MC results to illustrate 
the tightness of the deduced formulas. 

In Fig. |2l simulation results on the expectation of the 
MMSE, over 1.000 channel instances of the composite and 
the rain faded multibeam channels, are presented. The circle 
markers represent the simulated expectation over the channel 
instances, of the proposed approximation, given by (|9). On top 
of them, the analytic formulas presented in ( fT8b and (TToT l. are 
plotted. It is clear that these expressions precisely describe the 
expected values of the proposed approximation, thus providing 
a strong analytic tool for the description of the performance 
of the investigated systems. In the same figure, the tightness 
of the proposed approximation with respect to the actual 
performance, as calculated by MC simulations on (O, is also 
evident. Over the whole SNR region, the maximum deviation 
from the actual MMSE value is no more than 1.5dB for the 
composite fading case, while for the rain fading, the deviation 
is less than ldB. Especially for the SNR regions around 12.5 
dB and 2dB for the two cases respectively approximation 
becomes almost exact. In practical systems, these regions are 
of main interest. Consequently, the proposed expressions are 
very tight for systems with finite users and conventional link 
budgets, in contrast to asymptotic results based on large system 
dimensions or high SNR approaches. 



V. Conclusions 

A new, accurate approximation for the MMSE of a SIMO 
MU system that exhibits full receive correlation has been 
derived in the present paper, under an analytic framework that 
generalizes to various fading distributions. In particular, we 
have focused on the mobile composite satellite return link, 
as well as fixed satellite systems impaired by rain fading. We 
derive an analytic expression which approximates the expected 
value of the system performance in terms of MMSE, for finite 
system dimensions and over a large SNR region. 
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